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I describe the footprints of the classical chaos of the Yang-Mills fields in the quantum 
description. I also review the behavior of the BKL chaotic approach to the classical 
singularity on the basis of the Loop Quantum Gravity. 
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1. Introduction: Classical chaos in the Yang-Mills field theory and 
r*~^ General Relativity 

At the MG IX I had a talk on the chaos of the non-abelian gauge theory and 
the gravity (mainly cosmology). This talk was devoted to the classical fields. They 
Oh. are the gauge theories and arc non- linear. Thus, chaoticity is not surprising gen- 

(— i | erally,although we know examples of the stable solutions (solitons and others) to 

the non-linear equations. Concerning the Yang-Mills (YM) theory it is known that 
this theory is not integrable. All attempts to prove its integrability, not given up 
until now, are not conclusive and if one analysis them he will be convinced that 
the proof contains some conjecture (mostly of the mathematical nature). Of course 
some approximation scheme (e.g. large A-expansion) leads to the integrability. 

Therefore, it is not surprising that in 1981 it was found that the classical source- 
less YM equations exhibit the strong chaotic phenomena. This fact was proved ini- 
tially for the simplest form of the spatially homogeneous YM equations for n = 2, 3 
numbers of degrees of freedom [1,2] and then was extended to the spherically sym- 
metric field theory [3] . The spatial-temporal chaos of YM fields was established by 
the lattice calculations (see book [4] for review of this activity) . This gave rise to a 
program to describe the YM dynamics not in terms of the potential and fields but 
rather in terms of the loops (strings) variables [5]. 

These results demonstrate that the classical YM fields lack any special stable 
configuration, all states are chaotic and no particular configuration dominates in 
the Minkowski space-time in contrast to the Euclidean case, where the instantons 
give a dominant contribution to the functional integral. 

Turning now to the problem of the chaos in the General Relativity (GR) we have 
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to remark that despite an even longer history of the study the chaos in GR, the 
chaoticity is observed mostly in the problem of the approaching to the space like 
singularity t = 0. These studies are initiated and developing now due to the famous 
BKL chapter of GR [6]. The problem of the chaos in GR closely related to the 
problem of the singularity in the classical GR, where it is unavoidable according to 
the fundamental singularity theorems [7] . The backward evolution of an expanding 
universe leads to a singular state where the classical theory fails to be applied. 

If, as remarked by Landau (see [8]) one considers the metric g as a function of 
the synchronous time t only, at some finite time interval g = det(g flu ) tends to zero 
as t —* 0, independent of the equation of state or the character of the gravitational 
field, and this results the singularity. 

The corresponding metric ds 2 = dt 2 — dl 2 (a, (3 = 1,2,3), dl 2 — r ) a pdx a dx f} near 
the singularity t — is 

dl 2 = t 2pi dx\ + t 2p2 dx 2 2 + t 2p3 dx 2 3 (1) 

with 

i=i i=i 

This is so-called Kasner solution [9] corresponding to the Bianchi I spatial geom- 
etry. Thus one arrives to the regular flat, homogeneous, anisotropic space with the 
total volume homogeneously approaching the singularity. BKL consider the Bianchi 
IX geometry which generalize (1) and corresponds to the diagonal homogeneous 
anisotropic spatial metric with -fjj = a 2 

dl 2 = ajn{n I p dx a dx f3 (2) 

with the unit vectors along the axes n J a , (I = 1,2,3). a/ positive scale factors are 
functions of t only. 

The diagonal form of the matrix is a result of the vacuum Einstein equation 
Ro a = . However, even in the presence of matter, as argued by BKL, the pos- 
sible non-diagonal terms do not affect the character of the Kasner epochs and the 
character of the "replacement" of the Kasner exponents. As a result, the evolu- 
tion towards singularity proceeds via a series of successive oscillations during which 
the distances along two of the principal axes oscillate while they shrink monotoni- 
cally along the third axis (Kasner epochs). The volume V — J Vdet"fdx\dx2dx 3 — 
16ir 2 aiaiiani oc t as t — > 0. 

A new " era" begins when the monotonically falling metric components begin to 
oscillate while one of the previously oscillating directions begin to contract. This 
approach to the singularity reveals itself as an infinite succession of alternating 
Kasner "epochs" and has the character of a random process [10]. Thus infinite 
number of the oscillations are confined between any finite time t and moment of 
t = 0. 
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The central role in the BKL approach to the singularity plays the justification 
that the time derivatives dominate over space ones at the approach to the singular- 
ity. This fact allows to think that the inhomogeneous model can be well described 
by Bianchi IX model where the spatial geometry can be viewed as an assembly 
of the small patches each of which evolves almost independently. In other words, 
the dynamical decoupling of the different spatial points on the space-like slices has 
place. Sizes of the patches are defined by the scale of the space derivatives during 
the evolution while the curvatures grow. To sustain the homogeneity of the evolu- 
tion, parches have to be subdivide more and more at the vicinity of the singularity. 
If the geometry by some reason is discrete, such fragmentation must stop [11]. 

In the recent paper [12] this view was justified by the generalization the Misner's 
Mixmaster model to the generic inhomogeneous case. It is shown that neglect of 
spatial gradients possible in the asymptotic regime. In other words, authors of [12] 
claim that the generic cosmological solution near the singularity is isomorphic, point 
by point in space, to the one of the Bianchi VIII and IX models because the spatial 
coordinates in the Mixmaster model enter as parameters, and one can apply the 
long wavelength approximation. 

The global chaoticity of the classical YM fields and the chaotic behavior in the 
approach to the singularity in the classical GR, of course, have to be changed when 
the quantum effects enter the game. 

For the YM fields we have to take into account the quantum effects since QCD 
describes the quantum world of the hadrons and their interactions. For distances 
close to the Planck scale the gravitational field acquires large curvature and the 
evolution to the t = singularity has to be replaced by the quantum dynamics. 

In this situation, we expect that the chaos of the classical fields should be di- 
minished if not eliminated completely due to the quantum fluctuations of the gauge 
and quark fields. Below we show how it is happened and to which extent. We will 
see in general that the chaotic phenomena of the YM fields do not disappear in full, 
exhibiting the explicit footprints of the classical chaos in the quantum world. 



2. Quantum chaos of the YM fields 

The quantum insight into the YM dynamics clearly is achieved if we consider the 
spatially homogeneous potentials A®(t) (a = 1, 2, 3; = 0, 1, 2, 3) (YM classical and 
quantum mechanics). 

It is obvious that if this type of fields exhibits classical chaos, as we know, the 
wider class of the non- homogeneous YM fields A®(x, t) also will have this property. 
We mentioned that in the Introduction. 

The YM sourceless equation for potential A^(t) in the gauge A^t) = is reduced 
to the discrete Hamiltonian system (see [4] for review) 



d 2 A a Jdt 2 + g 2 {A a ] A b J A b l - A)A)A$) 



(3) 
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With the conserved "external" and "internal" angular momenta 



Mi = e llk A<*A a k , 



(4) 



N a = 



ahc A\Al, 



(5) 



which are vanished for the sourceless fields. Thus there exist seven integrals of 
motion and system is not integrable. 

Homogeneous limit of YM equations, or their long wavelength regime corre- 
sponds to the gluon high density n g and/or strong coupling regime g 2 n g \ >> 1. 
In this sense it is stated that the homogeneous fields are the relevant degrees of 
freedom for infrared regime [13]. Very recently, the equations (3) were obtained as 
a strong coupling limit of YM fields [14]. It was shown that in the leading order of 
1/g YM equations are reduced to (3). The resulting theory is stable and leads to 
the mass gap with the confinement. Author even calculated in this approach the 
glueball spectrum which is in good agreement with the lattice QCD computations. 

Consider now the simplest case of two and three degrees of freedom n = 2, 3 with 
x = A\,y = A\,z — A\. We obtain the systems of two (three) coupled oscillators 
with the potential energies 



Classically these systems, despite their extremely simple forms, exhibit strong 
chaotic behavior. Potential (6) {x 2 y 2 -model) has been used in various fields of 
science, including chemistry, astronomy, astrophysics and cosmology (chaotic in- 
flation). We mainly describe here the case n = 2. 

Quantum mechanical system with the potential (6) has only discrete spectrum 
[15] despite its open hyperbolic channels along the axes and its infinite phase space. 
Physically it is clear why it is so: quantum fluctuations, e.g. zero modes forbid the 
trajectory to escape along the axis where the potential energy vanishes. The system 
is thus confined to a finite volume and this implies the discreteness of the energy 
levels. Classically of course the "particle" always can escape along one of the axes 
without increasing its energy. 

Despite so drastic influence of the quantization on the behavior of the system 
(6) chaos left its footprints: periodic (unstable) orbits of the classical potential (6) 
after quantization have so called scars [16] (See also [17,18]). Energy level spacing 
distribution for the system (6) has the Wigncr-Dyson type distribution in contrast 
to the Poissonian one for the systems whose classical counterparts are regular. They 
are in accordance with the Random Matrix Theory for GOE. 

One can go further and propose that the traces of the classical chaos, in principle, 
should show in the real spectra of hadrons. For instance, if we would collect the 
rich enough glueball spectra then their mass spacing distribution has to reflect 
the chaoticity of the classical gluon field (gluon statistics has to be deal with the 



V(x,y) = (g 2 /2)x 2 y 2 



(6) 



V{x, y, z) = (g 2 /2) (x 2 y 2 + y 2 z 2 + x 2 z 2 ). 



(7) 
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assembly of the particles with the same quantum numbers) [19]. Not having today 
(when?) such a rich collection of glueballs as a cleanest sample, author of [20] used 
the relatively rich baryon and meson spectra for the examine the nearest-neighbor 
level spacing distributions for mass m < 2.5 GeV. It is seen that these distributions 
are well described by the Wigner surmise corresponding to the statistics of the GOE. 

Of course, one should consider this result as very preliminary since barions and 
mesons with their quark content are not the best case to check this idea (see [4] to 
inquire why the glueballs are necessary to solve this problem) . 

We would like to stress the important role which play here the so called billiards 
(classical and quantum). If we write the potential V(x 2 y 2 ) in the form {x 2 y 2 ) 1 ^ a 
where < a < 1, then limit a — corresponds to the so-called hyperbolic billiards, 
where the classical trajectory undergoes elastic collisions on an infinite barrier (hy- 
perbolic cylinder x 2 y 2 = 1). Trajectories lie in the x — y plane and consist of 
rectilinear segments constructed by the rules of geometrical optics. The notion of 
billiards plays an important role also in the GR [21]. 

Consider now the quantum mechanical adventures of the coupled YM oscillators 
in the study of the partition function for our non-integrable system (6) 

oo 

Z(t) = T r exp (-fcff)] = e ~ tE " ( 8 ) 



71=0 



with the quantum Hamiltonian H 



ry h 2 ( d 2 d 2 \ g 2 2 2 



2 \dx 2 dy 2 J 2 

using the quite effective method of the adiabatic separation of the motion in the 
hyperbola channels of the cquipotcntial curves xy = const [22]. 

The partition function defines the integrated density states N(E) by the inverse 
Laplace transform of Z(t) 

N(E) = dE p(E) = L- 1 (^j^j (10) 

and for the large enough energy levels E is given by the Thomas-Fermi term - the 
zero order term of the Wigner-Kirkwood expansion [23] 



Z (t) = _ 1 —— - (in + 9ln2 - C 



with C the Euler constant. 

From (11) and (8) one obtains N(E) <~ E 3 ^ 2 logE. For the hyperbola bil- 
liard (a = 0) the computations give [24], with the logarithmic precision, N(E) — 
±ElogE. 

This result some time ago was encouraging from the point of view of the famous 
Hilbert-Polya-Berry program to look for the quantum (classically chaotic) Hamil- 
tonian whose eigenvalues reproduce the Ricmann zeta-function's zeros. However, 
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from Random Matrix Theory we know that such Hamiltonian (or some operator ) 
must be non-invariant with respect to time inversion (GUE not GOE!). 

Returning to the calculation of Z(t), we apply the method of [22] (see also [25]): 
The range of the integration over x and y variables and momenta p x and p y in 
the calculation of Z(t) in the Winger representation is divided into two regions: 
the central region (Jx/Jy/ < Q , where Q is arbitrary, not specified scale) and 
the channels region (Q < /x/, Q < /y/).In the central region it is natural to apply 
Wigncr-Kirkwood method; in the channels with the "slow" motion in the x variables 
(along the channel) again will be used the Wigner-Kirkwood expansion, but the 
"fast" motion, transverse to the channel, must be treated quantum mechanically. 
Remarkably, the dependence on Q from both regions is cancelled in the leading 
terms (tQ 4 )^ 1 << 1 up to eighth order of h. No doubts are left that the higher 
order of h terms behave similarly (see [25]). Thus, only Q-indcpendent, non- leading 
terms are contributing to Z(t). 

We encounter here the phenomenon of the "transmutation" of small parameters: 
Classical parameter 1/tQ 4 which rules the adiabatic separation of the variables in 
the channels transmutes into the small quantum parameter of the final asymptotic 
series for Z(t). 

It is interesting to follow the motion in the channel (along the x axis) a little in 
detail. Motion in the channel can be describe by the Hamiltonian 



where w = gx is x dependent frequency and eigenvalues of (12) are (n + 1/2) hgx. 

In the channel (/x/ >> /y/) where the derivatives w .r. t. x are small relative 
to the derivatives w. r .t. y, we may first average the motion over the quantum 
fluctuations of y [25,26] described by (12) and by the corresponding wave func- 
tion involving Hermite polynomials with the frequency uj — gx. The corresponding 
average 



then appears as an effective potential for the motion in the "slow" variable x 



This is the well known Schrodingcr equation for a linear potential having the 
solutions in terms of Airy functions which shows the linear confinement and discrete 
spectrum for eigenvalues. 

We would like to emphasize that this confinement is not like standard phe- 
nomenon commonly refcred to as quark confinement. Here the potential described by 
the gauge field amplitude x(t). One may call this phenomenon as "self-confinement" : 
the fields themselves "prepare" the effective potential barrier, prohibiting escape to 
the infinity. 



(12) 



< n/Hy/n >= (n + -) hgx 




(13) 
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As we remarked above, just due to this consistent treatment of the motion 
in the channels when each n-th quantum evolves along the x axis according the 
Hamiltonian 

H { x n) = \P 2 X + (n+ 1/2) hg\x\ (14) 

there are the precise cancellation of all leading quantum corrections (in the regime 
<< 1) and only non-leading but Q-independent corrections survive and lead 
to the final answer for Z(t) in the form of asymptotic series with the expansion 
parameter g 2 h 4 t 2 . 

One remark is worthy on the discreteness of the spectrum. In the corresponding 
supersymmetric quantum mechanics, due to the cancellation between the bosonic 
and fermionic modes, there appears the continuous spectrum which coexists with 
the discrete one [27] and the confinement generally has not place. 

There is another approach to calculate Z(t): add to the potential (6) the extra 
term V 2 (x 2 + y 2 ) - Higgs vacuum term, compute Z and then put V — [28]. This 
limit opens the hyperbola channels classically. In quantum mechanics, this limit 
leads to the singularities: Logarithmic for the Thomas-Fermi term and power like 
for the higher quantum corrections V for the fc-th order of h. However, quantum 
mechanics cures that: it introduce the Higgs-like term [29] to the potential and due 
to this the limit V = in Z(T) has no singularity at all. 

In conclusion of this Section we state that the lessons derived from this quantum 
mechanical study of the higher order corrections to the homogeneous limit of YM 
equations would be useful for the better understanding the internal dynamics of the 
YM quantum field theory. 

3. BKL chaos in the Loop Quantum Gravity (Loop Quantum 
Cosmology) 

In this Section we consider which kind modifications on should expect for the os- 
cillating chaotic approach to the classical singularity t — (BKL scenario) if one 
includes the quantum corrections. It is natural that at each novel scheme of the 
quantum effects to the gravitational field, BKL chapter is tested with the various 
conclusion and verdicts about this scenario. I described this efforts briefly in the 
Talk to MG IX mentioned in the Introduction [30] . Now the number of these consid- 
erations of the BKL behavior is increased significantly and includes Matrix models, 
String theory, its brane aspects, higher derivative corrections, matter content modi- 
fications (dilaton, p-form fields) etc. (see e.g., [31]). In the most of these approaches 
classical singularity, as a rule, are not avoided and this has a strong influence on 
the BKL scheme. 

Here we describe only one approach based on the Loop Quantum Gravity (LQG) 
[32] and its sibling, Loop Quantum Cosmology (LQC) which is based not only on 
the LQC, but includes several additional assumptions and simplifications [33]. 

Here we describe very briefly the basic notions of LQG and LQC. 
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As it is known, this approach to the canonical Hamiltonian G R is based on the 
use not the ADM variables (spatial metric and extrinsic curvature) but the Ashtckar 
variables [34,35] based on the inclusion the spin connection variables what allows 
the formulation closely to the YM-like gauge field theory. 

Briefly, the following has a place. In the expression of the contra variant spatial 
metric g ab = ef e\ in terms of triad (orthogonal and normalized at each point) there 
is redundancy due to an arbitrary 3-dimensional rotation of the triad which does 
not change metric. 

Densitizcd triad 

Et=g 1/2 e1, (g = detg ab ) (15) 

together with the SU(2) connection A l a (Ashtekar connection), j = 1, 2, 3; a, b = 
1, 2, 3) form the pair (A l a , Ef) canonically connected to the metric conjugated pair 
(g a b, p ab )- As in the gauge theories, "momentum" Ef of the spin connection A l a 
is an analog of the electric field of SU(2) YM theory (index labels a basic element 
of the SU(2) Lie algebra). Connections A % a involve the curvature of space and spin 
connection T l a 

Ai = r a +l3Ki (16) 

where K l a ("torsion") defines the extrinsic curvature K ai = K a i,e b . Positive param- 
eter (3 was introduced by Barbero [36] as substitute of the former imaginary unit, 
to have a certain reality conditions for Ashtekar variables. 

This /3-ambiguity leads to the Poisson bracket dependent on [3: 

{A l a {x),E b J (y)} = nP5 h J)5{x,y) (k = 8ttG) (17) 

and after the Dirac first class constraint quantization, leads to the /3-dependent 
Hamiltonian constraint which rules the evolution of the system ADM-like way. 

The rest two constraint: Gauss constraint (generating triad rotations) and dif- 
feomorphism constraint (generating spatial diffeomorphisms) are independent on 
the Barbero-Immirzi ambiguity. 

We would like to remark that parameter (3 can be considered as the rescaling of 
the triad 

K = -^9< (18) 

and this allows the possibility to introduce the conformal symmetry [37]. 

However, in the Dirac quantization this will lead to the new first class constraint 
corresponding to the conformal symmetry. 

The main advantage of the new variables is that they allow a natural smearing 
of the basis fields (A, E) to the linear objects without introducing a background and 
retain the well-defined algebra. The connections integrated along a curve, exponen- 
tiated with a path-ordered way. Thus we arrive to the holonomies (this is analogous 
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of the quantum mechanics where the Heisenberg operator, e.g., x is represented by 
Weyl operator e lx ): 



Oi are the Pauli matrices. 

Trace of h e [A] corresponds to the Wilson loop for the closed curve in the YM 
theory. Similarly, we arrive to the fluxes by integrating densitized triad over two- 
surfaces S: 



where n a is the normal to the surface S. 

The above introduced smearing, without introducing background, eliminate all 
delta functions in the Poisson relations giving the well-defined algebra to construct 
the Hilbert space. We do not dwell on the fundamental problem of LQG of this 
construction in terms of the holonomies and fluxes in the conditions of the diffco- 
morphism invariancc. There are wide spectrum of the conceptual and the technical 
problems for LQG (see e.g. [38]). 

We only make some remarks: 

* Barbieri- Immirzi ambiguity can be resolved " experimentally" comparing 
the LQG results with the Bekenstein-Hawking entropy for the Black Hole 
[39] . Taking the formula of the area eigenvalues A — 8ir[3lp \J j(j + 1) one 
obtains (3 = ln2/irV3. The new data yield 0.27398. 

* ft was remarked [40] that the role of parameter (3 in the canonical quantum 
gravity is analogous in various senses to that of the parameter describing 
the different sectors associated with the topological structure of the finite 
gauge transformations in the YM theory. In contrast to the LQG, the 0-term 
enters as the Pontryagin topological term which is a total derivative. 

* After quantization, homologies and fluxes act as well defined operators, 
fluxes have the discrete spectra. Since the spatial geometry is determined by 
the densitized triads, spatial geometry is discrete as well, with the discrete 
area and volume operators. In this sense, it is sometime declared that the 
Quantum Gravity is a Natural Regulator of matter [41]. By this reason, 
differential equation for the Wheeler-De Witt constraint is replaced by the 
difference equations. 

* Although in the LQG there is some parallel with the Wilson loops of 
YM theory (hence the term Loop Quantum Gravity), there is essential 
difference. In the YM theory different sizes of loops are inequivalent in the 
light of the interpretation relying to the quark confinement (e.g., inside the 
large loops fields are chaotic, inside the small loops are regular, small-large 
w.r.t . the confinement radius). The value of the Wilson loop is invariant 




(19) 



1 




(20) 
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under continuum deformations only for the vanishing field strength. In the 
LQG, with its diffcomorphism invariance, there is no physical information 
in the shape and size since two networks of the different shape but the same 
topology can always be related by a suitable diffcomorphism, independently 
of the "value" of the Ashtekar field strength. 
After these remarks we turn to the BKL problem for Bianchi IX space geometry 

We will base here on the mini-superspace spanned by the scale factors ai of (2), 

diagonal anisotropic model, gn = aj (t). 

Hamiltonian constraint is written in terms of the spin connections T/ (below we 

follow [42]). 



(21) 



2 I 

H = - [(TjT K - Ti) aj - - aiajiiK + {TiTj - T K ) a K 
k 4 

~ a K aiaj + {TkTt - Tj) a.j - i ajd^d/] 

1 faj a K aj \ 1 fp K p J p J p K \ 

2 \a K + a r aja K J 2 (p^ J 

(I, J, K) an even permutation of 1, 2, 3; p 1 = e IKL ax o-l ■ 

Derive now the classical equations transforming to a new canonical variables 717 
and q 1 and to new time coordinate: 

717 = -(log ai) , Q ,/ = ^logp / , dt = a 1 a 2 a 3 dT 

2 (23) 

W, Tj} = K 5 j , 

Separating terms with momenta like variables 717, we obtain the potential term 



W(ai, 02, 03) = ^ 



— a\ a\ — a\ a\ — a\a\. (24) 



From (22) and (23) it is seen that at the small ai (or p 1 ) due to the divergencies 
of the spin connection components there is singularity. 
Classical equation of motion are 

^(logai)" = {al-4f-a\ (25) 

and two e.o.m. by cyclic. 

Right hand sides in (25) are vanished for Bianchi I geometry, giving a/ ~ t ai 
with, = I = aj, i.e. the Kasner solution. 

For Bianchi IX geometry one may write the evolution potential in terms oip 1 : 

w( P \ V \ p 3 ) - 2 [(pV (r! r 2 - r 3 ) +p i P 1 (r^r, - r 2 ) 

(26) 

+p 2 p 3 (r 2 r 3 -ri)] 
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which has an infinite walls at p 1 w due to the divergence of the spin connec- 
tion components. Evolution consists of the succession of the Kasner epochs with 
reflections on the walls and this process never stops what leads to the BKL chaos. 
To be closer to the common picture one can introduce the Misner variables 

1 1 

^ = -3 lo S^ = -3 log(aia 2 a 3 ) 

and anisotropies /?/: 

Then the potential (26) takes the form 

W(fl, p+, /?_) = \e- in [e- sl3 + -4e- 2 ?+ cosh(2\/3 /3_) 

2 (27) 

2e 4/3 +(cosh(4\/3/3_)-l)]. 

Volume dependence factorizes and the rest anisotropy potential shows exponen- 
tial walls for the large anisotropies. For instance, at typical wall has a form if one 
takes (3_ = and (3 + < 0: 

W~^e- in ~ s P+ (28) 

To prevent the "eternal" reflections at the walls where the expansion/contraction 
behavior of the different directions change, it is necessary to stop the unconstrained 
rise of the heights of the walls. Just the quantum effects are called up for this 
prevention. 

In other words, they should lead to the upper limit on the curvature. 

What is the concrete scenario to achieve this aim in the LQC? 

First of all, one has to have in mind that in the LQG and in the LQC there exist 
effective minimal length (or area, Ai/ 2 — Stt (3l 2 y/3/2) or the maximal curvature. 

Quantization according to the rules of the game in LQG results in the replace- 
ment of the spin connections components Ti by the effective coefficients, which leads 
to the effective potential instead of (26). 

Central moment here is the special rules of the quantization of the inverse den- 
sitized triad variables (p 1 ) -1 or the inverse volume, giving that in the LQG they 
are not singular at p 1 = despite the classical curvature divergence. 

One should remember that the LQC which actively considers various important 
phcnomcnological effects is not in the strong sense the direct limiting case of the full 
LQG where the desired boundedness of the inverse scale factors or the inverse vol- 
ume are ensured. LQC is the usual cosmological mini superspace with its symmetry 
reduction and oversimplifications, quantized by the LQG methods and techniques. 
For this reason, for instance isotropic model in the LQC, in contrast to some claims, 
has no bounded from above inverse scale factor whereas in the full scale LQG it is 
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proved that such a inverse scale and inverse volume operators have bounded from 
above eigenvalues [43] . The physical explanation of this not common situation may 
be that the isotropic homogeneous quantum fluctuations for isotropic mini super- 
space model are not enough to eliminate the classical singularity. For anisotropic 
Bianchi IX model it is not excluded that the fluctuations of the same symmetry 
may ensure this elimination although it is not based on the firm grounds as it has 
place for the LQG [43]. 

Anyway, taking the assumption that for Bianchi IX geometry this conjecture 
is realized, let us continue to follow what is happened with the chaos near the 
singularity. We again follow [42] . For convenience, we take | (3 l 2 = 1 making p 1 
dimensionless. 

Quantization replaces {p l )~ l in the spin connection components by the function 
F(p/2j) where the parameter j appears explicitly and controls the peak of the 
function F. The same parameter j enters the expression of the area in the LQG: 
A(j) = 8 ir (3 lp^/j (j + 1). Further, follow the recipe of the LQG to obtain the 
effective description, one needs to replace all negative powers of triad p 1 with the 
appropriate factors of the spectrum of the inverse volume operator. 

For instance, 

p-^^d = D(p/p*)/ P 3/2 

with p* = a* 2 = 16 7r j fi n /3. Function D s=a 1 for p/p* >> 1, recovering the classical 
behavior. For the small p (or a), D{^-) ~ p x l 2 or d ~ p 6 thus giving the smooth 
behavior at the singularity. 

For the anisotropic homogeneous model the components of (p 1 )^ 1 (I = 1,2,3) 
are replaced in the spin connection components by a function F{p ! /2 j) giving the 
effective spin connection. Parameter j belongs to the set of the ambiguities of this 
framework and controls the peak of the function F. The resulting effective potential 
as a function of p 1 at fixed volume V has a form 



A 



i 



p 'P > [pi 



= W j (jp\p\ 2jq)^ 

(29) 



V 4 F 2 (n) 

w [3 - 2qF(q)] 



where q= ± 



At the peak and beyond it F(q) ~l/q and we have the classical wall \ e ^ 4 °^ 8 ' 3 + . 
The peak of the finite walls is reached for a constant q which in the usual variables 
gives that e - 2n + 2 P+ = const. Maxima of the wall lie on the line (3+ — ft + const in 
the classical phase space and the height of the wall drops off as e~ mi <~ V 4 with 
the decreasing V — > 0. 

At very small volume the walls collapse more rapidly and the effective potential 
becomes negative everywhere at the volumes close of in the Planck units. For the 
smallest value of j = 1/2 it is about Planck volume. 
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Thus, with the decreasing walls during the evolution towards the singularity the 
classical reflections will stop at a finite time interval and the chaos should disappear. 
Universe -at some time- can "jump over the wall" [42], Kasner regime becomes 
stable. If we think about non homogeneities, the patches of the corresponding space 
become of the order of a Planck volume, i.e. the scale of the discreteness. Below 
that scale further fragmentation does not happen and the discreteness is preserved. 

In the large volumes when the evolution is chaotic, two nearby points - patches 
will diverge away with no correlations between the points ("non interacting two 
dimensional gas" ) . In the vicinity of singularity, the chaotic motion is replaced by the 
Kasner evolution, points-patches begin to correlate ("interacting two dimensional 
gas"). 

The evolution to the singularity on the basis of LQC with its final non-chaotic 
scenario near the Planck scale and beyond rises the important question of the in- 
creasing the role of non-homogeneities at that scale. To sustain the homogeneous 
regime one needs the further and further fragmentation of patches of the spatial 
regions. But at the conditions of the discrete spatial geometry, the fragmentation 
must be stopped and the homogeneity should be replaced by non - homogeneity. 
Thus, at the approaching to the singularity role of the inhomogeneous quantum 
fluctuations may be essential. This brings us to the relatively old notion of the 
"turbulent" universe [50], [51]. 

4. In lieu of conclusion 

In lieu of conclusion, we enumerate here several important problems considered by 
the LQC and not concerning the chaos in the cosmology. LQC, using the similar 
approach (based on the minisuperspace and the recipes of the quantization from the 
full LQG) has a several important contribution to the " explanation" of the inflation 
[44] , to the quantum nature of the Big Bang [45] , possible observational signatures 
in the CMBR [46], avoidance of the future singularity [47] where the interesting 
effect of the negative quadratic density correction inspired by LQC is observed in 
the FRW equation, and the quantum evaporation of the naked singularity [48] . Last 
paper gives an interesting view on the problem of the gravitational collapse of the 
matter (scalar field as an example) near the classical singularity. The authors of [48] 
observed the rise of the strong outward energy flux which dissolves the collapsing 
cloud before the formation of the singularity. This effect based on the LQC may be 
considered as a mechanism of the censorship of the naked singularity [49]. 

Authors of [48] think about the observational signature of this effect in the 
astrophysical bursts. Time will show how reliable are these interesting investigations. 
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